Sub-femtotesla scalar atomic magnetometer using multipass cells 



(N 

o 

(N 

< 



Of 



o 

-I— > . 
cd. 

Of 



>: 

On, 

OV 

O' 
t-h ; 

00' 

o: 



X 



D. Sheng, 1 S. Li, 1 ' 2 N. Dural, 1 and M. V. Romans 1 

1 Department of physics, Princeton University, Princeton, New Jersey 08544, USA 
2 Department of Measurement Technology and Instruments, 
Zhejiang University of Science and Technology, Hangzhou 310023, China 
(Dated: August 7, 2012) 

Scalar atomic magnetometers have many attractive features but their sensitivity has been relatively poor. 
We describe a Rb scalar gradiometer using two multi-pass optical cells. We use a pump-probe measure- 
ment scheme to suppress spin-exchange relaxation and two probe pulses to find the spin precession zero 
crossing times with a resolution of 1 psec. We realize magnetic field sensitivity of 0.54 fT/Hz 1 / 2 , which 
improves by an order of magnitude the best scalar magnetometer sensitivity and surpasses the quantum 
limit set by spin-exchange collisions for a scalar magnetometer with the same measurement volume oper- 
ating in a continuous regime. 
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Alkali-metal magnetometers can surpass SQUIDs as the 
most sensitive detectors of magnetic field, reaching sensi- 
tivity below 1 fT/Hz 1 / 2 H H, but only if they are oper- 
ated near zero magnetic field to eliminate spin relaxation 
due to spin-exchange collisions Many magnetome- 

ter applications, such as searches for permanent electric 
dipole moments [5], detection of NMR signals 0], and 
low-field magnetic resonance imaging [0], require sensi- 
tive magnetic measurements in a finite magnetic field. In 
addition, scalar magnetometers measuring the Zeeman fre- 
quency are unique among magnetic sensors in being in- 
sensitive to the direction of the field, making them par- 
ticularly suitable for geomagnetic mapping 18|] and field 
measurements in space H [loll . The sensitivity of scalar 
magnetometers has been relatively poor, as summarized re- 
cently in II 1111 . The best directly measured scalar magne- 
tometer sensitivity is equal 7 fT/Hz 1//2 with a measurement 
volume of 1.5 cm 3 1121 . while estimates of fundamental 
sensitivity per unit measurement volume for various types 
of scalar alkali-metal magnetometers range from several 
fT crn^/Hz 1 / 2 mB to about 1 fT crn^/Hz 1 / 2 HI. 



Here we describe a new type of scalar atomic magnetome- 
ter using multi-pass vapor cells lfl5L [Irjll and operating in 
a pulsed pump-probe mode [fTvh to achieve magnetic field 
sensitivity of 0.54 fT/Hz 1 / 2 with a measurement volume of 
0.66 cm 3 in each multi-pass cell. The magnetometer sensi- 
tivity approaches, for the first time, the fundamental limit 
set by Rb-Rb collisions. We also develop here a quantita- 
tive method to analyze significant effects of atomic diffu- 
sion on the spectrum of the spin-projection noise in vapor 
cells with buffer gas using a spin time-correlation function. 

The sensitivity of an atomic magnetometer, as any other 
frequency measurement, is fundamentally limited by spin 
projection noise and spin relaxation [18]. For N spin- 
1/2 atoms with coherence time T 2 the sensitivity after 
a long measurement time f > T 2 is given by 5B = 
y/2e/NT 2 t/'y, where 7 is the gyromagnetic ratio. Spin 
squeezing techniques can reduce this uncertainty by a fac- 



tor of yfe, but do not change the scaling with iV 118I42C 
The number of atoms can be increased until collisions be- 
tween them start to limit T 2 . Writing T 2 X = nav, where 
n is the density of atoms, a is the spin relaxation cross- 
section, and v is the average collisional velocity, and taking 
t = 0.5 sec to calculate the magnetic field spectral noise 
density B n in T/Hz 1 / 2 , we obtain 



B r , 



(2/7)^^^. 
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Thus the magnetic field spectral noise density per mea- 
surement volume V is fundamentally limited by the spin- 
relaxation cross-section. It also sets the limit on the mini- 
mum energy resolution per unit bandwidth e = -B 2 V /2^ 
of atomic magnetometers, which can, in certain cases, ap- 
proach H [2J]. In hot alkali-metal vapor magnetometers op- 
erating in a finite magnetic field the relaxation is dominated 
by the spin-exchange cross-section a$E = 1-9 x 10~ 14 
cm 3 , which gives from Eq. (TJ a limit of about 1.3 fT 
crn^/Hz 1 / 2 for 7 = 2vr x 700 kHz/G. 

However, alkali-metal spin-exchange is a nonlinear pro- 
cess, which can modify the fundamental sensitivity given 
by Eq. (H). The magnetic resonance linewidth can be nar- 
rowed by optical pumping of atoms into a stretched spin 
state Bill , but fundamental sensitivity for a scalar magne- 
tometer still remains limited by the spin-exchange cross- 
section if it is operated in a continuous optical pumping 
regime Ill2ll . The limit calculated in [12] taking into ac- 
count 87 Rb nuclear spin I = 3/2 is 0.51 fT crn^/Hz 1 / 2 . 
On the other hand, if the magnetometer is operated in a 
pulsed pump-probe regime and takes advantage of spin- 
squeezing techniques, the sensitivity can be asymptotically 
limited by the spin-destruction cross-section, which is as 
low as a SD = 10~ 18 cm 3 for K atoms, leading to a poten- 
tial improvement by 2 orders of magnitude Il22ll . Thus, it 
is particularly interesting to study spin projection noise in 
scalar alkali-metal magnetometers, both because it presents 
a real limit to their practical sensitivity and because of large 
improvement possible from spin-squeezing techniques. 
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FIG. 1. (Color online) (a) Experiment setup. PBS: polarized 
beam splitter, PMF: polarization maintaining fiber, DAQ: data ac- 
quisition card, (b) The timing of the pulsed operation, (c) Optical 
rotation (black line) recorded for one probe pulse at atom density 
of 0.8 x 10 13 /cm 3 together with a fitted curve (red dash line), (d) 
Magnetic field noise spectrum obtained in the gradiometer mode 
in presence of a calibrating gradient magnetic field at 40 Hz. 



A key parameter for measurements of spin-projection 
noise is the optical depth on resonance OD = a nl, where 
(To i s probe laser absorption cross-section on resonance and 
I is the path length of the probe beam through the atomic 
vapor [23]. We have developed multi-pass optical cells 
with mirrors internal to the alkali-metal vapor cell to in- 
crease I by two orders of magnitude [16]. Compared to 
optical cavities, multi-pass cells have a much larger inter- 
action volume and allow direct recording of large optical 
rotations. We use two 42-pass cells placed in the same va- 
por cell as a gradiometer with a baseline equal to the 1 .5 
cm distance between the cells, see Fig.QIa). The cells have 
cylindrical mirrors with 10 cm radius of curvature sepa- 
rated by 30 mm. One of the mirrors in each cell has a 2.5 
mm diameter hole for entrance and exit of the probe beam 
focused to a waist diameter of 1.9 mm. The glass vapor 
cell contains a drop of enriched 87 Rb and 70 torr N 2 gas. 
A boron-nitride oven is used to heat the vapor cell using 
AC currents at 600 kHz to 120°C, giving an OD ~ 5000. 
The cell is placed in a bias magnetic field of 72.9 mG in 
the z direction generated by an ultra-stable custom current 
source and is enclosed in a 5-layer magnetic shield. 

We measure the atom density n from the transverse re- 
laxation T 2 at low polarization, which is dominated by 
spin-exchange collisions with a known cross-section lfl6ll . 
The number of atoms participating in the measurement at 
any given time N = nVf, is determined from the area 
of Faraday rotation power spectral density for unpolar- 
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FIG. 2. (a) Unpolarized atom spin noise spectrum taken at 0.42 
G with atom density of 1.2xl0 13 /cm 3 after subtraction of the 
photon shot noise background from a spin noise measurement 
at low magnetic field, (b) Spin correlation function due to the 
diffusion of the atoms out of the probe beam. Solid line (dashed 
line) is the experimental (calculated) result with atom density of 
1.2x 10 13 /cm 3 . The inset shows the calculated beam pattern at 
the center of the cavity. 

ized atoms Il24ll . We make measurements of the noise 
peak at two different magnetic fields and take their differ- 
ence to remove the background dominated by photon shot 
noise. Fig. |2a) shows one example of unpolarized power 
spectral density obtained using this method, which gives 
V b = 0.35(2) cm 3 for each cell. 

While diffusion does not affect the area under the spin 
noise peak, it causes the lineshape of the noise spectrum to 
deviate from a simple Lorentzian. To analyze it quantita- 
tively, we consider the time autocovariance function of the 
Faraday rotation signal (f)(t), which is given by the Fourier 
transform of the power spectrum. One can show that 

n (Ff) J /( ri )G( ri - r 2 , r)/(r 2 )d 3 r l( Z 3 r 2 , (2) 

where the sum is taken over the two alkali-metal hyper- 
fine states, F a = I + 1/2 and F b = I - 1/2, and 
{Ff) = F i (F l + l)(2F + 1)/6(2J + 1). The dis- 
persion factor Di = l/(^, — v) for far detuning of the 
probe frequency v from the hyperfine resonances Z/,-. Here 
I(r) is the total probe laser intensity at position r, includ- 
ing all beam passes inside the cavity, and G(r, r) is the 
Green's function for spin evolution with a diffusion coef- 
ficient D and a transverse relaxation time T 2 , G(r, r) = 

e -r 2 /iDr-r /T 2 J ^ Dr y/2 _ j ntensity profiIe of me 

probe laser in the cell is determined by measuring the input 



Gaussian beam size and calculating the astigmatic Gaus- 
sian beam propagation in the multi-pass cell 12511 . An ex- 
ample of the calculated intensity profile in the middle of the 
cell is shown in the inset of Fig. 2. The effective number 
of atoms participating in the measurement is defined as the 
number of atoms that would generate the same spin noise 
area (</>(i) 2 ) if interrogated with a uniform probe intensity. 
We obtain a generalization of a result given in Ref . lE4ll that 
works for laser beams with varying focusing and overlap, 

Ul{r)dyd Z f 

n /J(r)W ' {i) 

where / is the total probe laser path in the multi-pass 
cell. Based on calculated intensity profile we obtain 
Vfc = 0.36 cm 3 , in good agreement with direct experi- 
mental measurements. In Fig. 2 we compare the diffusion 
component of the calculated spin time-correlation function 
C d ( T ) = (0(t)</>(t + r))e t / T V(0(t) 2 ) with the experi- 
mental measurement obtained by taking the Fourier trans- 
form of the spin noise peak after centering it at zero fre- 
quency and correcting for the transverse spin relaxation 
time C,j(t) = C(r)e r / T2 . They agree well except at early 
times due to deviations from a perfect Gaussian of tightly 
focused beams within the cell, indicating that Green's func- 
tion method can quantitatively describe the lineshape of the 
spin noise spectrum in the presence of diffusion with mul- 
tiple overlapping laser beams. 

FigureQJb) shows the timing for magnetic field measure- 
ments. For optical pumping pulse, which lasts 14 msec, 
we use two circular polarized beams on resonance with the 
Dl transitions from both ground hyperfine states. Then 
we apply a tt/2 rf pulse lasting 3 periods of the Zeeman 
resonance frequency. We apply the first probe light pulse 
shortly after the rf excitation and the second probe pulse 
with a delay time T from the first one. The probe laser 
is tuned to 794.780 nm and the power of the light exiting 
from each multi-pass cell is about 0.5 mW. We turn on and 
off the probe light slowly compared with the Larmor pe- 
riod using an AOM to suppress transient spin excitation. 
The pump-probe cycle is repeated every 16.6 msec, syn- 
chronized with 60 Hz to reduce its influence. 

Figure [Tic) shows a typical record of the optical rotation 
signal during one of the probe pulses. We fit the data using 
the equation 11611 : 

1' - .Isin ( 2<> ( 1 - t -^ L ) sin(w(t - t c )) + ip )+B. 



(4) 

We find the time of zero crossings t c i,t c2 of the first and 
second pulse and calculate T c = t c2 — t c i which gives 
a measure of the magnetic field B = 2Trm/jT c , where 
m is the number of spin precession cycles between mea- 
sured zero crossings. If the measurements are repeated 
with overall duty cycle d, then the magnetic field sensi- 
tivity per Hz 1 / 2 is given by: 
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B n = B5T c J2/dT c , 



(5) 



FIG. 3. The dependence of the standard deviation 5T C on the rf 
pulse amplitude. Empty (solid) circles denote the back-action 
evasion with stroboscopic modulation (normal modulation) of 
probe light. Here the probe pulse length and T are equal to four 
Larmor periods. 



where 5T C is the standard deviation of repeated measure- 
ments of T c . 

The two multi-pass cells work as a gradiometer to mea- 
sure dB z I dy with a noise level which is \[2 larger than 
given by Eq. (0. While a scalar magnetometer does not 
require calibration, we check its response by applying a 
known source of magnetic field gradient. Fig. [He) shows 
Fourier transform of repeated magnetic field measurements 
in the presence of a gradient dB z /dy with rms ampli- 
tude of 21.6 fT/cm oscillating at 40 Hz. It introduces 
a 33 fT magnetic field difference between the centers of 
the two cells. For this measurement the atomic density 
is 1.4x 10 13 /cm 3 , with the probe pulse length of four Lar- 
mor periods, the separation between two probe pulses T of 
823 [is, and the cycle period is 5 ms. The integration of 
40 Hz peak in the spectrum gives the rms signal of 32.4 fT 
and confirms the sensitivity of the gradiometer. 

The limiting fundamental noise sources include spin pro- 
jection noise (SPN) and photon shot noise (PSN), while 
technical sources include magnetic shield noise and time 
jitter of the data acquisition. One important feature of 
our arrangement is back-action evasion of quantum fluc- 
tuations of the probe beam circular polarization due to zero 
spin polarization of atoms in the z direction following the 
tt/2 pulse. Fig. [3] shows the dependence of the noise on 
the rf excitation amplitude when it deviates from the tt/2 
amplitude. We compare it to the noise in T c when us- 
ing a stroboscopic probe modulation back-action evasion 
scheme 122J], where the probe beam is modulated at twice 
the Larmor frequency with a duty cycle 20% in the probe 
pulse period. The results confirm that the magnetometer 
works in back-action-free regime. The magnetic shield gra- 
dient noise is due by thermal Johnson currents and is calcu- 
lated based on known electrical conductivity of the shields, 
giving 0.40(5) fT/Hz 1 / 2 [26]. The time jitter noise is deter- 
mined by recording the signals from the same multi-pass 
cell with two acquisition channels and ranges from 0.3 to 
0.5 ps depending on the length of the probe pulse. 
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Fig. [4] shows the magnetic field sensitivity for a range 
of parameters. The experimental duty cycle is limited to 
10% by available pump laser power, but we use d = 1 in 
Eq. (O to find the fundamental sensitivity from the mea- 
sured uncertainty 5T C . Fig.Ua) shows the nonlinear relax- 
ation of transverse spin polarization due to spin exchange 
at four different densities, from which we find that the ini- 
tial transverse polarization is equal to 0.96(1). We plot 
the sensitivity as a function of the probe pulse length t p 
in Fig. IHb). The variance in T c due to PSN and data ac- 
quisition noise decreases as l/t p and the variance due to 
SPN also decreases because atom diffusion effectively in- 
volves more atoms into the measurement. The effective 
number of atoms N m participating in the measurement af- 
ter a pulse time t p can be found using the diffusion correla- 
tion function N m = nV b / 2/t p j* p (l - t/t p )C d (t)dt . 

For the longest pulse length of 230 /isec we obtain V m = 
1.9Vf,, corresponding to an effective interaction volume 
of 0.66 cm 3 . We also show theoretical estimate of the 
sensitivity including ASN, PSN, magnetic gradient noise 
and time jitter noise in Fig. Ub) with solid lines and only 
ASN and PSN with broken lines. Figure Uc) shows sim- 
ilar results at other densities. When the atom density in- 
creases, the optimal T decreases because of faster spin 
relaxation, indicating that the magnetometer works in Rb 
collision-limited regime. For the longest probe pulse length 
and atom density of 1.4x 10 13 /cm 3 , the experimental data 
shows a best sensitivity of 0.54 fT/Hz 1 / 2 , which is 10% 
above the predicted value. In the absence of magnetic 
shield noise the fundamental sensitivity is projected to be 
0.3 fT/Hz 1 / 2 , dominated by ASN. 

In conclusion, we described a scalar magnetometer 
based on multipass atomic vapor cells. It uses a pulsed 
mode with a high initial polarization and reaches the spin- 
exchange collision limited regime where the sensitivity is 
largely independent of atom density. The best sensitivity 
obtained is 0.54 fT/Hz 1 / 2 with an effective interaction vol- 
ume of 0.66 cm 3 , which is an order of magnitude improve- 
ment over the previous best sensitivity for a scalar magne- 
tometer and exceeds the quantum-limited sensitivity of a 
scalar magnetometer operating in a continuous mode. We 
also developed a quantitative method for analyzing the ef- 
fect of diffusion on quantum spin noise using spin time- 
correlation function. By relying on precision timing mea- 
surements with very wide dynamic range and fractional 
field sensitivity of 7 x 10 _11 /Hz 1 / 2 this magnetometer 
opens the possibility of fundamentally new applications, 
for example unshielded detection of magnetoencephalog- 
raphy signals II27I1 . The sensitivity per unit volume can 
be further improved in this system by reducing the decay 
of the spin time-correlation function due to atomic diffu- 
sion, which will allow suppression of ASN due to spin- 
squeezing between two probe pulses. The spin correlation 
decay is dominated by a few tightly focused beam spots in 
the multi-pass cell and can be reduced by modifying multi- 
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FIG. 4. (color) (a) Relaxation of transverse polarization Pt 
(points) for n =1.9 (green), 1.4 (blue), 0.8 (red) and 0.6 (black) 
x 10 13 /cm 3 together with theoretical prediction (lines), (b) The 
points are the experimental results of magnetic field sensitivity 
at n=1.4x 10 13 /cm 3 , with probe pulse length of four (orange), 
eight (cyan) and twelve (purple) Larmor periods. Solid (dashed) 
lines are the theoretically calculated sensitivity based on mea- 
sured parameters with (without) magnetic shield noise and time 
jitter noise, (c) Experiment results (points) of magnetic field sen- 
sitivity with probe length of twelve Larmor periods at different 
atom densities, with the same color notation as plot (a), together 
with theoretical predictions (lines). 



pass cell parameters to avoid tight beam focusing. This 
work was supported by DARPA. 
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